SherbyNet - /page/wikibook/entry.mu

Micron Document

<!DOCTYPE html>

<title>Data assimilation</title>

</head>

Data assimilation

</h1>

</header>

/* start https://en.wikipedia.org/ */

.mw-parser-output .ambox{border:1px solid #a2a9b1;border-left:10px solid #36c;background-color:#fbfbfb;box-sizing:border-box}.mw-parser-output .ambox+link+.ambox,.mw-parser-output .ambox+link+style+.ambox,.mw-parser-output .ambox+link+link+.ambox,.mw-parser-output .ambox+.mw-empty-elt+link+.ambox,.mw-parser-output .ambox+.mw-empty-elt+link+style+.ambox,.mw-parser-output .ambox+.mw-empty-elt+link+link+.ambox{margin-top:-1px}html body.mediawiki .mw-parser-output .ambox.mbox-small-left{margin:4px 1em 4px 0;overflow:hidden;width:238px;border-collapse:collapse;font-size:88%;line-height:1.25em}.mw-parser-output .ambox-speedy{border-left:10px solid #b32424;background-color:#fee7e6}.mw-parser-output .ambox-delete{border-left:10px solid #b32424}.mw-parser-output .ambox-content{border-left:10px solid #f28500}.mw-parser-output .ambox-style{border-left:10px solid #fc3}.mw-parser-output .ambox-move{border-left:10px solid #9932cc}.mw-parser-output .ambox-protection{border-left:10px solid #a2a9b1}.mw-parser-output .ambox .mbox-text{border:none;padding:0.25em 0.5em;width:100%}.mw-parser-output .ambox .mbox-image{border:none;padding:2px 0 2px 0.5em;text-align:center}.mw-parser-output .ambox .mbox-imageright{border:none;padding:2px 0.5em 2px 0;text-align:center}.mw-parser-output .ambox .mbox-empty-cell{border:none;padding:0;width:1px}.mw-parser-output .ambox .mbox-image-div{width:52px}@media(min-width:720px){.mw-parser-output .ambox{margin:0 10%}}@media print{body.ns-0 .mw-parser-output .ambox{display:none!important}}

/* end https://en.wikipedia.org/ */

</style><style data-mw-deduplicate="TemplateStyles:r1248332772">

/* start https://en.wikipedia.org/ */

.mw-parser-output .multiple-issues-text{width:95%;margin:0.2em 0}.mw-parser-output .multiple-issues-text>.mw-collapsible-content{margin-top:0.3em}.mw-parser-output .compact-ambox .ambox{border:none;border-collapse:collapse;background-color:transparent;margin:0 0 0 1.6em!important;padding:0!important;width:auto;display:block}body.mediawiki .mw-parser-output .compact-ambox .ambox.mbox-small-left{font-size:100%;width:auto;margin:0}.mw-parser-output .compact-ambox .ambox .mbox-text{padding:0!important;margin:0!important}.mw-parser-output .compact-ambox .ambox .mbox-text-span{display:list-item;line-height:1.5em;list-style-type:disc}body.skin-minerva .mw-parser-output .multiple-issues-text>.mw-collapsible-toggle,.mw-parser-output .compact-ambox .ambox .mbox-image,.mw-parser-output .compact-ambox .ambox .mbox-imageright,.mw-parser-output .compact-ambox .ambox .mbox-empty-cell,.mw-parser-output .compact-ambox .hide-when-compact{display:none}

/* end https://en.wikipedia.org/ */

</style>

Data assimilation refers to a large group of methods that update information from numerical computer models with information from observations. Data assimilation is used to update model states, model trajectories over time, model parameters, and combinations thereof. What distinguishes data assimilation from other estimation methods is that the computer model is a dynamical model, i.e. the model describes how model variables change over time, and its firm mathematical foundation in Bayesian Inference. As such, it generalizes inverse methods and has close connections with machine learning.

Data assimilation initially developed in the field of <a href="Numerical_weather_prediction" title="Numerical weather prediction">numerical weather prediction</a>. Numerical weather prediction models are equations describing the evolution of the atmosphere, typically coded into a computer program. When these models are used for forecasting the model output quickly deviates from the real atmosphere. Hence, we use observations of the atmosphere to keep the model on track. Data assimilation provides a very large number of practical ways to bring these observations into the models.

Simply inserting point-wise measurements into the numerical models did not provide a satisfactory solution. Real world measurements contain errors both due to the quality of the instrument and how accurately the position of the measurement is known. These errors can cause instabilities in the models that eliminate any level of skill in a forecast. Thus, more sophisticated methods were needed in order to initialize a model using all available data while making sure to maintain stability in the numerical model. Such data typically includes the measurements as well as a previous forecast valid at the same time the measurements are made. If applied iteratively, this process begins to accumulate information from past observations into all subsequent forecasts.

Because data assimilation developed out of the field of numerical weather prediction, it initially gained popularity amongst the geosciences. In fact, one of the most cited publication in all of the geosciences is an application of data assimilation to reconstruct the observed history of the atmosphere.<a href="#cite_note-1">[1]</a>

<div class="mw-heading mw-heading2"><h2 id="Details_of_the_data_assimilation_process">Details of the data assimilation process</h2></div>

Classically, data assimilation has been applied to chaotic dynamical systems that are too difficult to predict using simple extrapolation methods. The cause of this difficulty is that small changes in initial conditions can lead to large changes in prediction accuracy. This is sometimes known as the <a href="Butterfly_effect" title="Butterfly effect">butterfly effect</a> – the sensitive dependence on <a href="Initial_condition" title="Initial condition">initial conditions</a> in which a small change in one state of a <a href="Deterministic" class="mw-redirect" title="Deterministic">deterministic</a> <a href="Nonlinear_system" title="Nonlinear system">nonlinear system</a> can result in large differences in a later state.

At any update time, data assimilation usually takes a forecast (also known as the first guess, or background information) and applies a correction to the forecast based on a set of observed data and estimated errors that are present in both the observations and the forecast itself. The difference between the forecast and the observations at that time is called the departure or the <a href="Innovation_(signal_processing)" title="Innovation (signal processing)">innovation</a> (as it provides new information to the data assimilation process). A weighting factor is applied to the innovation to determine how much of a correction should be made to the forecast based on the new information from the observations. The best estimate of the state of the system based on the correction to the forecast determined by a weighting factor times the innovation is called the analysis. In one dimension, computing the analysis could be as simple as forming a weighted average of a forecasted and observed value. In multiple dimensions the problem becomes more difficult. Much of the work in data assimilation is focused on adequately estimating the appropriate weighting factor based on intricate knowledge of the errors in the system.

The measurements are usually made of a real-world system, rather than of the model's incomplete representation of that system, and so a special function called the observation operator (usually depicted by h() for a nonlinear operator or "H" for its linearization) is needed to map the modeled variable to a form that can be directly compared with the observation.

<div class="mw-heading mw-heading2"><h2 id="Data_assimilation_as_statistical_estimation">Data assimilation as statistical estimation</h2></div>

One of the common mathematical philosophical perspectives is to view data assimilation as a Bayesian estimation problem. From this perspective, the analysis step is an application of <a href="Bayes'_theorem" title="Bayes' theorem">Bayes' theorem</a> and the overall assimilation procedure is an example of <a href="Recursive_Bayesian_estimation" title="Recursive Bayesian estimation">recursive Bayesian estimation</a>. However, the probabilistic analysis is usually simplified to a computationally feasible form. Advancing the probability distribution in time would be done exactly in the general case by the <a href="Fokker%E2%80%93Planck_equation" title="Fokker–Planck equation">Fokker–Planck equation</a>, but that is not feasible for high-dimensional systems; so, various approximations operating on simplified <a href="Representation_(mathematics)" title="Representation (mathematics)">representations</a> of the probability distributions are used instead. Often the probability distributions are assumed <a href="Normal_distribution" title="Normal distribution">Gaussian</a> so that they can be represented by their mean and covariance, which gives rise to the <a href="Kalman_filter" title="Kalman filter">Kalman filter</a>.

Many methods represent the probability distributions only by the mean and input some pre-calculated covariance. An example of a direct (or sequential) method to compute this is called optimal statistical interpolation, or simply optimal interpolation (OI). An alternative approach is to iteratively solve a cost function that solves an identical problem. These are called "variational methods", such as 3D-Var and 4D-Var. Typical minimization algorithms are the <a href="Conjugate_gradient_method" title="Conjugate gradient method">conjugate gradient method</a> or the <a href="Generalized_minimal_residual_method" title="Generalized minimal residual method">generalized minimal residual method</a>. The <a href="Ensemble_Kalman_filter" title="Ensemble Kalman filter">ensemble Kalman filter</a> is sequential method that uses a Monte Carlo approach to estimate both the mean and the covariance of a Gaussian probability distribution by an ensemble of simulations. More recently, hybrid combinations of ensemble approaches and variational methods have become more popular (e.g. they are used for operational forecasts both at the <a href="European_Centre_for_Medium-Range_Weather_Forecasts" title="European Centre for Medium-Range Weather Forecasts">European Centre for Medium-Range Weather Forecasts</a> (ECMWF) and at the <a href="NOAA" class="mw-redirect" title="NOAA">NOAA</a> <a href="National_Centers_for_Environmental_Prediction" title="National Centers for Environmental Prediction">National Centers for Environmental Prediction</a> (NCEP).

<div class="mw-heading mw-heading2"><h2 id="Data_assimilation_as_a_model_update">Data assimilation as a model update</h2></div>

Data assimilation can also be achieved within a model update loop, where we will iterate an initial model (or initial guess) in an optimisation loop to constrain the model to the observed data. Many optimisation approaches exist and all of them can be set up to update the model, for instance, <a href="Evolutionary_algorithm" title="Evolutionary algorithm">evolutionary algorithm</a> have proven to be efficient as free of hypothesis, but computationally expensive.

<div class="mw-heading mw-heading2"><h2 id="Weather_forecasting_applications">Weather forecasting applications</h2></div>

In numerical weather prediction applications, data assimilation is most widely known as a method for combining observations of meteorological variables such as <a href="Temperature" title="Temperature">temperature</a> and <a href="Atmospheric_pressure" title="Atmospheric pressure">atmospheric pressure</a> with prior forecasts in order to initialize numerical forecast models.

<div class="mw-heading mw-heading3"><h3 id="Necessity">Necessity</h3></div>

The <a href="Atmosphere" title="Atmosphere">atmosphere</a> is a <a href="Fluid" title="Fluid">fluid</a>. The idea of numerical weather prediction is to sample the state of the fluid at a given time and use the equations of <a href="Fluid_dynamics" title="Fluid dynamics">fluid dynamics</a> and <a href="Thermodynamics" title="Thermodynamics">thermodynamics</a> to estimate the state of the fluid at some time in the future. The process of entering observation data into the model to generate <a href="Initial_value_problem" title="Initial value problem">initial conditions</a> is called initialization. On land, terrain maps available at resolutions down to 1 kilometer (0.6 mi) globally are used to help model atmospheric circulations within regions of rugged topography, in order to better depict features such as downslope winds, <a href="Lee_wave" title="Lee wave">mountain waves</a> and related cloudiness that affects incoming solar radiation.<a href="#cite_note-2">[2]</a> The main inputs from country-based weather services are observations from devices (called <a href="Radiosonde" title="Radiosonde">radiosondes</a>) in weather balloons that measure various atmospheric parameters and transmits them to a fixed receiver, as well as from <a href="Weather_satellite" title="Weather satellite">weather satellites</a>. The <a href="World_Meteorological_Organization" title="World Meteorological Organization">World Meteorological Organization</a> acts to standardize the instrumentation, observing practices and timing of these observations worldwide. Stations either report hourly in <a href="METAR" title="METAR">METAR</a> reports,<a href="#cite_note-3">[3]</a> or every six hours in <a href="SYNOP" title="SYNOP">SYNOP</a> reports.<a href="#cite_note-4">[4]</a> These observations are irregularly spaced, so they are processed by data assimilation and objective analysis methods, which perform quality control and obtain values at locations usable by the model's mathematical algorithms.<a href="#cite_note-Krishnamurti_Annu_Rev_FM-5">[5]</a> Some global models use <a href="Finite_differences" class="mw-redirect" title="Finite differences">finite differences</a>, in which the world is represented as discrete points on a regularly spaced grid of latitude and longitude;<a href="#cite_note-6">[6]</a> other models use <a href="Spectral_methods" class="mw-redirect" title="Spectral methods">spectral methods</a> that solve for a range of wavelengths. The data are then used in the model as the starting point for a forecast.<a href="#cite_note-7">[7]</a>

A variety of methods are used to gather observational data for use in numerical models. Sites launch radiosondes in weather balloons which rise through the <a href="Troposphere" title="Troposphere">troposphere</a> and well into the <a href="Stratosphere" title="Stratosphere">stratosphere</a>.<a href="#cite_note-8">[8]</a> Information from weather satellites is used where traditional data sources are not available. Commerce provides <a href="Pilot_report" title="Pilot report">pilot reports</a> along aircraft routes<a href="#cite_note-9">[9]</a> and ship reports along shipping routes.<a href="#cite_note-10">[10]</a> Research projects use <a href="Weather_reconnaissance" title="Weather reconnaissance">reconnaissance aircraft</a> to fly in and around weather systems of interest, such as <a href="Tropical_cyclone" title="Tropical cyclone">tropical cyclones</a>.<a href="#cite_note-Hurricane_Hunters-11">[11]</a><a href="#cite_note-SunHerald-12">[12]</a> Reconnaissance aircraft are also flown over the open oceans during the cold season into systems which cause significant uncertainty in forecast guidance, or are expected to be of high impact from three to seven days into the future over the downstream continent.<a href="#cite_note-13">[13]</a> Sea ice began to be initialized in forecast models in 1971.<a href="#cite_note-14">[14]</a> Efforts to involve <a href="Sea_surface_temperature" title="Sea surface temperature">sea surface temperature</a> in model initialization began in 1972 due to its role in modulating weather in higher latitudes of the Pacific.<a href="#cite_note-15">[15]</a>

<div class="mw-heading mw-heading3"><h3 id="History">History</h3></div>

In 1922, <a href="Lewis_Fry_Richardson" title="Lewis Fry Richardson">Lewis Fry Richardson</a> published the first attempt at forecasting the weather numerically. Using a <a href="Hydrostatic_equilibrium" title="Hydrostatic equilibrium">hydrostatic</a> variation of <a href="Vilhelm_Bjerknes" title="Vilhelm Bjerknes">Bjerknes's</a> <a href="Primitive_equations" title="Primitive equations">primitive equations</a>,<a href="#cite_note-Lynch_JCP-16">[16]</a> Richardson produced by hand a 6-hour forecast for the state of the atmosphere over two points in central Europe, taking at least six weeks to do so.<a href="#cite_note-Lynch_Ch1-17">[17]</a> His forecast calculated that the change in <a href="Atmospheric_pressure" title="Atmospheric pressure">surface pressure</a> would be 145 <a href="Bar_(unit)" title="Bar (unit)">millibars</a> (4.3 <a href="Inch_of_mercury" title="Inch of mercury">inHg</a>), an unrealistic value incorrect by two orders of magnitude. The large error was caused by an imbalance in the pressure and wind velocity fields used as the initial conditions in his analysis,<a href="#cite_note-Lynch_JCP-16">[16]</a> indicating the need for a data assimilation scheme.

Originally "subjective analysis" had been used in which numerical weather prediction (NWP) forecasts had been adjusted by meteorologists using their operational expertise. Then "objective analysis" (e.g. Cressman algorithm) was introduced for automated data assimilation. These objective methods used simple interpolation approaches, and thus were 3DDA (three-dimensional data assimilation) methods.

Later, 4DDA (four-dimensional data assimilation) methods, called "nudging", were developed, such as in the <a href="MM5_(weather_model)" title="MM5 (weather model)">MM5</a> model. They are based on the simple idea of Newtonian relaxation (the 2nd axiom of Newton). They introduce into the right part of dynamical equations of the model a term that is proportional to the difference of the calculated meteorological variable and the observed value. This term that has a negative sign keeps the calculated <a href="State_space_(controls)" class="mw-redirect" title="State space (controls)">state vector</a> closer to the observations. Nudging can be interpreted as a variant of the <a href="Kalman-Bucy_filter" class="mw-redirect" title="Kalman-Bucy filter">Kalman-Bucy filter</a> (a continuous time version of the <a href="Kalman_filter" title="Kalman filter">Kalman filter</a>) with the gain matrix prescribed rather than obtained from covariances.

A major development was achieved by L. Gandin (1963) who introduced the "statistical interpolation" (or "optimal interpolation") method, which developed earlier ideas of Kolmogorov. This is a 3DDA method and is a type of <a href="Regression_analysis" title="Regression analysis">regression analysis</a> which utilizes information about the spatial distributions of <a href="Covariance" title="Covariance">covariance</a> functions of the errors of the "first guess" field (previous forecast) and "true field". These functions are never known. However, the different approximations were assumed.

The optimal interpolation algorithm is the reduced version of the <a href="Kalman_filtering" class="mw-redirect" title="Kalman filtering">Kalman filtering</a> (KF) algorithm and in which the covariance matrices are not calculated from the dynamical equations but are pre-determined in advance.

Attempts to introduce the KF algorithms as a 4DDA tool for NWP models came later. However, this was (and remains) a difficult task because the full version requires solution of the enormous number of additional equations (~N*N~10**12, where N=Nx*Ny*Nz is the size of the state vector, Nx~100, Ny~100, Nz~100 – the dimensions of the computational grid). To overcome this difficulty, approximate or suboptimal Kalman filters were developed. These include the <a href="Ensemble_Kalman_filter" title="Ensemble Kalman filter">Ensemble Kalman filter</a> and the Reduced-Rank Kalman filters (RRSQRT).<a href="#cite_note-18">[18]</a>

Another significant advance in the development of the 4DDA methods was utilizing the <a href="Optimal_control" title="Optimal control">optimal control</a> theory (variational approach) in the works of Le Dimet and Talagrand (1986), based on the previous works of J.-L. Lions and G. Marchuk, the latter being the first to apply that theory in the environmental modeling. The significant advantage of the variational approaches is that the meteorological fields satisfy the dynamical equations of the NWP model and at the same time they minimize the functional, characterizing their difference from observations. Thus, the problem of constrained minimization is solved. The 3DDA variational methods were developed for the first time by Sasaki (1958).

As was shown by Lorenc (1986), all the above-mentioned 4DDA methods are in some limit equivalent, i.e. under some assumptions they minimize the same <a href="Loss_function" title="Loss function">cost function</a>. However, in practical applications these assumptions are never fulfilled, the different methods perform differently and generally it is not clear what approach (Kalman filtering or variational) is better. The fundamental questions also arise in application of the advanced DA techniques such as convergence of the computational method to the global minimum of the functional to be minimised. For instance, cost function or the set in which the solution is sought can be not convex. The 4DDA method which is currently most successful<a href="#cite_note-19">[19]</a><a href="#cite_note-20">[20]</a> is hybrid incremental 4D-Var, where an ensemble is used to augment the climatological background error covariances at the start of the data assimilation time window, but the background error covariances are evolved during the time window by a simplified version of the NWP forecast model. This data assimilation method is used operationally at forecast centres such as the <a href="Met_Office" title="Met Office">Met Office</a>.<a href="#cite_note-21">[21]</a><a href="#cite_note-22">[22]</a>

<div class="mw-heading mw-heading3"><h3 id="Cost_function">Cost function</h3></div>

The process of creating the analysis in data assimilation often involves minimization of a <a href="Loss_function" title="Loss function">cost function</a>. A typical cost function would be the sum of the squared deviations of the analysis values from the observations weighted by the accuracy of the observations, plus the sum of the squared deviations of the forecast fields and the analyzed fields weighted by the accuracy of the forecast. This has the effect of making sure that the analysis does not drift too far away from observations and forecasts that are known to usually be reliable.

</mrow>

</mrow>

<msub>

</mrow>

</mrow>

</msub>

<msup>

</mrow>

</msup>

<msup>

</mrow>

</mrow>

</msup>

</mrow>

<msub>

</mrow>

</mrow>

</msub>

</mrow>

</mrow>

</mrow>

<msup>

</mrow>

</msup>

<msup>

</mrow>

</mrow>

</msup>

</mrow>

</mrow>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle J(\mathbf {x} )=(\mathbf {x} -\mathbf {x} _{b})^{\mathrm {T} }\mathbf {B} ^{-1}(\mathbf {x} -\mathbf {x} _{b})+(\mathbf {y} -{\mathit {H}}[\mathbf {x} ])^{\mathrm {T} }\mathbf {R} ^{-1}(\mathbf {y} -{\mathit {H}}[\mathbf {x} ]),}</annotation>

</semantics>

</math><img src="./c5e0355c8f4568b6a4b73afb07eb8d9c53e704cc.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.838ex; width:60.488ex; height:3.176ex;" alt="{\displaystyle J(\mathbf {x} )=(\mathbf {x} -\mathbf {x} _{b})^{\mathrm {T} }\mathbf {B} ^{-1}(\mathbf {x} -\mathbf {x} _{b})+(\mathbf {y} -{\mathit {H}}[\mathbf {x} ])^{\mathrm {T} }\mathbf {R} ^{-1}(\mathbf {y} -{\mathit {H}}[\mathbf {x} ]),}" loading="lazy">

where <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle \mathbf {B} }">

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \mathbf {B} }</annotation>

</semantics>

</math><img src="./cafb0ef39b0f5ffa23c170aa7f7b4e718327c4d1.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.338ex; width:1.901ex; height:2.176ex;" alt="{\displaystyle \mathbf {B} }" loading="lazy"> denotes the background error covariance, <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle \mathbf {R} }">

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \mathbf {R} }</annotation>

</semantics>

</math><img src="./5de85fcc2a00d8ba14aae84aeef812d7fef4b3d5.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.338ex; width:2.003ex; height:2.176ex;" alt="{\displaystyle \mathbf {R} }" loading="lazy"> the observational error covariance.

</mrow>

<msup>

</mrow>

</mrow>

</msup>

</mrow>

<msub>

</mrow>

</mrow>

</msub>

<msup>

</mrow>

</mrow>

</msup>

<msup>

</mrow>

</mrow>

</msup>

</mrow>

</mrow>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \nabla J(\mathbf {x} )=2\mathbf {B} ^{-1}(\mathbf {x} -\mathbf {x} _{b})-2{\mathit {H}}^{T}\mathbf {R} ^{-1}(\mathbf {y} -{\mathit {H}}[\mathbf {x} ])}</annotation>

</semantics>

</math><img src="./e88d5d464b06b0488c4af88cd420a8dfefdd7d0c.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.838ex; width:45.479ex; height:3.176ex;" alt="{\displaystyle \nabla J(\mathbf {x} )=2\mathbf {B} ^{-1}(\mathbf {x} -\mathbf {x} _{b})-2{\mathit {H}}^{T}\mathbf {R} ^{-1}(\mathbf {y} -{\mathit {H}}[\mathbf {x} ])}" loading="lazy">

</mrow>

</mrow>

<msub>

</mrow>

</mrow>

</msub>

<msup>

</mrow>

</msup>

<msup>

</mrow>

</mrow>

</msup>

</mrow>

<msub>

</mrow>

</mrow>

</msub>

</mrow>

</mrow>

</munderover>

<msub>

</mrow>

</mrow>

</msub>

<msub>

</mrow>

</mrow>

</msub>

<msub>

</mrow>

</mrow>

</msub>

<msup>

</mrow>

</msup>

</mrow>

</mrow>

</mrow>

</msubsup>

<msub>

</mrow>

</mrow>

</msub>

<msub>

</mrow>

</mrow>

</msub>

<msub>

</mrow>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle J(\mathbf {x} )=(\mathbf {x} -\mathbf {x} _{b})^{\mathrm {T} }\mathbf {B} ^{-1}(\mathbf {x} -\mathbf {x} _{b})+\sum _{i=0}^{n}(\mathbf {y} _{i}-{\mathit {H}}_{i}[\mathbf {x} _{i}])^{\mathrm {T} }\mathbf {R} _{i}^{-1}(\mathbf {y} _{i}-{\mathit {H}}_{i}[\mathbf {x} _{i}])}</annotation>

</semantics>

</math><img src="./31d47e6941f24f5f489b798049325e8206ffbc16.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -3.005ex; width:67.994ex; height:6.843ex;" alt="{\displaystyle J(\mathbf {x} )=(\mathbf {x} -\mathbf {x} _{b})^{\mathrm {T} }\mathbf {B} ^{-1}(\mathbf {x} -\mathbf {x} _{b})+\sum _{i=0}^{n}(\mathbf {y} _{i}-{\mathit {H}}_{i}[\mathbf {x} _{i}])^{\mathrm {T} }\mathbf {R} _{i}^{-1}(\mathbf {y} _{i}-{\mathit {H}}_{i}[\mathbf {x} _{i}])}" loading="lazy">

provided that <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle {\mathit {H}}}">

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle {\mathit {H}}}</annotation>

</semantics>

</math><img src="./8c1274a23bea898e8c8743252319bb9808edbdde.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.338ex; margin-right: -0.272ex; width:1.999ex; height:2.176ex;" alt="{\displaystyle {\mathit {H}}}" loading="lazy"> is a linear operator (matrix).

<div class="mw-heading mw-heading3"><h3 id="Future_development">Future development</h3></div>

Factors driving the rapid development of data assimilation methods for NWP models include:

<ul><li>Utilizing the observations currently offers promising improvement in <a href="Forecast_skill" title="Forecast skill">forecast skill</a> at a variety of spatial scales (from global to highly local) and time scales.</li>

<li>The number of different kinds of available observations (<a href="Sodar" title="Sodar">sodars</a>, <a href="Radar" title="Radar">radars</a>, <a href="Satellite" title="Satellite">satellite</a>) is rapidly growing.</li></ul>

<div class="mw-heading mw-heading2"><h2 id="Other_applications">Other applications</h2></div>

<div class="mw-heading mw-heading3"><h3 id="Monitoring_water_and_energy_transfers">Monitoring water and energy transfers</h3></div>

Data assimilation has been used, in the 1980s and 1990s, in several HAPEX (Hydrologic and Atmospheric Pilot Experiment) projects for monitoring energy transfers between the soil, vegetation and atmosphere. For instance:

- <a rel="nofollow" class="external text" href="https://www.eol.ucar.edu/field_projects/hapex-mobilhy">HAPEX-MobilHy</a>,<a href="#cite_note-24">[24]</a> HAPEX-Sahel,<a href="#cite_note-25">[25]</a>

- the "Alpilles-ReSeDA" (Remote Sensing Data Assimilation) experiment,<a href="#cite_note-26">[26]</a><a href="#cite_note-27">[27]</a> a European project in the <a rel="nofollow" class="external text" href="https://cordis.europa.eu/programme/rcn/471/en">FP4-ENV</a> program<a href="#cite_note-28">[28]</a> which took place in the <a href="Alpilles" title="Alpilles">Alpilles</a> region, South-East of France (1996–97). The Flow-chart diagram (right), excerpted from the final report of that project,<a href="#cite_note-ReSeDA_report-23">[23]</a> shows how to infer variables of interest such as canopy state, radiative fluxes, environmental budget, production in quantity and quality, from remote sensing data and ancillary information. In that diagram, the small blue-green arrows indicate the direct way the models actually run.<a href="#cite_note-29">[29]</a>

<div class="mw-heading mw-heading3"><h3 id="Other_forecasting_applications">Other forecasting applications</h3></div>

Data assimilation methods are currently also used in other environmental forecasting problems, e.g. in <a href="Hydrology" title="Hydrology">hydrological</a> and hydrogeological forecasting.<a href="#cite_note-30">[30]</a> Bayesian networks may also be used in a data assimilation approach to assess natural hazards such as landslides.<a href="#cite_note-31">[31]</a>

Given the abundance of spacecraft data for other planets in the <a href="Solar_System" title="Solar System">Solar System</a>, data assimilation is now also applied beyond the Earth to obtain re-analyses of the atmospheric state of extraterrestrial planets. Mars is the only extraterrestrial planet to which data assimilation has been applied so far. Available spacecraft data include, in particular, retrievals of temperature and dust/water/ice optical thicknesses from the <a href="Thermal_Emission_Spectrometer" title="Thermal Emission Spectrometer">Thermal Emission Spectrometer</a> onboard NASA's <a href="Mars_Global_Surveyor" title="Mars Global Surveyor">Mars Global Surveyor</a> and the Mars Climate Sounder onboard NASA's <a href="Mars_Reconnaissance_Orbiter" title="Mars Reconnaissance Orbiter">Mars Reconnaissance Orbiter</a>. Two methods of data assimilation have been applied to these datasets: an Analysis Correction scheme <a href="#cite_note-montabone-32">[32]</a> and two Ensemble Kalman Filter schemes,<a href="#cite_note-hoffman-33">[33]</a><a href="#cite_note-wrf-34">[34]</a> both using a global circulation model of the martian atmosphere as forward model. The Mars Analysis Correction Data Assimilation (MACDA) dataset is publicly available from the British Atmospheric Data Centre.<a href="#cite_note-badc-35">[35]</a>

Data assimilation is a part of the challenge for every forecasting problem.

Dealing with biased data is a serious challenge in data assimilation. Further development of methods to deal with biases will be of particular use. If there are several instruments observing the same variable then intercomparing them using <a href="Probability_density_function" title="Probability density function">probability distribution functions</a> can be instructive.

The numerical forecast models are becoming of higher resolution due to the increase of <a href="Computational_power" class="mw-redirect" title="Computational power">computational power</a>, with operational atmospheric models now running with horizontal resolutions of order of 1 km (e.g. at the German National Meteorological Service, the <a href="Deutscher_Wetterdienst" title="Deutscher Wetterdienst">Deutscher Wetterdienst</a> (<abbr>DWD</abbr>) and <a href="Met_Office" title="Met Office">Met Office</a> in the UK). This increase in horizontal resolutions is starting to allow to resolve more chaotic features of the non-linear models, e.g. to resolve <a href="Convection" title="Convection">convection</a> on the grid scale, or clouds, in the atmospheric models. This increasing non-linearity in the models and <a href="Inverse_problem" title="Inverse problem">observation operators</a> poses a new problem in the data assimilation. The existing data assimilation methods such as many variants of <a href="Ensemble_Kalman_filter" title="Ensemble Kalman filter">ensemble Kalman filters</a> and variational methods, well established with linear or near-linear models, are being assessed on non-linear models.

Many new methods are being developed, e.g. <a href="Particle_filter" title="Particle filter">particle filters</a> for high-dimensional problems, and hybrid data assimilation methods.<a href="#cite_note-36">[36]</a>

Other uses include trajectory estimation for the <a href="Apollo_program" title="Apollo program">Apollo program</a>, <a href="GPS" class="mw-redirect" title="GPS">GPS</a>, and <a href="Autochem" title="Autochem">atmospheric chemistry</a>.

<ul><li><a href="Calibration_(statistics)" title="Calibration (statistics)">Calibration</a></li></ul>

<div class="mw-heading mw-heading2"><h2 id="References">References</h2></div>

/* start https://en.wikipedia.org/ */

.mw-parser-output .reflist{margin-bottom:0.5em;list-style-type:decimal}@media screen{.mw-parser-output .reflist{font-size:90%}}.mw-parser-output .reflist .references{font-size:100%;margin-bottom:0;list-style-type:inherit}.mw-parser-output .reflist-columns-2{column-width:30em}.mw-parser-output .reflist-columns-3{column-width:25em}.mw-parser-output .reflist-columns{margin-top:0.3em}.mw-parser-output .reflist-columns ol{margin-top:0}.mw-parser-output .reflist-columns li{page-break-inside:avoid;break-inside:avoid-column}.mw-parser-output .reflist-upper-alpha{list-style-type:upper-alpha}.mw-parser-output .reflist-upper-roman{list-style-type:upper-roman}.mw-parser-output .reflist-lower-alpha{list-style-type:lower-alpha}.mw-parser-output .reflist-lower-greek{list-style-type:lower-greek}.mw-parser-output .reflist-lower-roman{list-style-type:lower-roman}

/* end https://en.wikipedia.org/ */

</style><div class="reflist reflist-columns references-column-width" style="column-width: 30em;">

/* start https://en.wikipedia.org/ */

.mw-parser-output cite.citation{font-style:inherit;word-wrap:break-word}.mw-parser-output .citation q{quotes:"\"""\"""'""'"}.mw-parser-output .citation:target{background-color:rgba(0,127,255,0.133)}.mw-parser-output .id-lock-free.id-lock-free a{background:url("./mw/Lock-green.svg")right 0.1em center/9px no-repeat}.mw-parser-output .id-lock-limited.id-lock-limited a,.mw-parser-output .id-lock-registration.id-lock-registration a{background:url("./mw/Lock-gray-alt-2.svg")right 0.1em center/9px no-repeat}.mw-parser-output .id-lock-subscription.id-lock-subscription a{background:url("./mw/Lock-red-alt-2.svg")right 0.1em center/9px no-repeat}.mw-parser-output .cs1-ws-icon a{background:url("./mw/Wikisource-logo.svg")right 0.1em center/12px no-repeat}body:not(.skin-timeless):not(.skin-minerva) .mw-parser-output .id-lock-free a,body:not(.skin-timeless):not(.skin-minerva) .mw-parser-output .id-lock-limited a,body:not(.skin-timeless):not(.skin-minerva) .mw-parser-output .id-lock-registration a,body:not(.skin-timeless):not(.skin-minerva) .mw-parser-output .id-lock-subscription a,body:not(.skin-timeless):not(.skin-minerva) .mw-parser-output .cs1-ws-icon a{background-size:contain;padding:0 1em 0 0}.mw-parser-output .cs1-code{color:inherit;background:inherit;border:none;padding:inherit}.mw-parser-output .cs1-hidden-error{display:none;color:var(--color-error,#d33)}.mw-parser-output .cs1-visible-error{color:var(--color-error,#d33)}.mw-parser-output .cs1-maint{display:none;color:#085;margin-left:0.3em}.mw-parser-output .cs1-kern-left{padding-left:0.2em}.mw-parser-output .cs1-kern-right{padding-right:0.2em}.mw-parser-output .citation .mw-selflink{font-weight:inherit}@media screen{.mw-parser-output .cs1-format{font-size:95%}html.skin-theme-clientpref-night .mw-parser-output .cs1-maint{color:#18911f}}@media screen and (prefers-color-scheme:dark){html.skin-theme-clientpref-os .mw-parser-output .cs1-maint{color:#18911f}}

/* end https://en.wikipedia.org/ */

</style><cite id="CITEREFKalnayand_coauthors1996" class="citation journal cs1">Kalnay, Eugenia; and coauthors (1996). <a rel="nofollow" class="external text" href="https://doi.org/10.1175%2F1520-0477%281996%29077%3C0437%3ATNYRP%3E2.0.CO%3B2">"The NCEP/NCAR 40-Year Reanalysis Project"</a>. Bulletin of the American Meteorological Society. 77 (March): 437–471. <a href="Bibcode_(identifier)" class="mw-redirect" title="Bibcode (identifier)">Bibcode</a>:<a rel="nofollow" class="external text" href="https://ui.adsabs.harvard.edu/abs/1996BAMS...77..437K">1996BAMS...77..437K</a>. <a href="Doi_(identifier)" class="mw-redirect" title="Doi (identifier)">doi</a>:<a rel="nofollow" class="external text" href="https://doi.org/10.1175%2F1520-0477%281996%29077%3C0437%3ATNYRP%3E2.0.CO%3B2">10.1175/1520-0477(1996)077<0437:TNYRP>2.0.CO;2</a>. <a href="ISSN_(identifier)" class="mw-redirect" title="ISSN (identifier)">ISSN</a> <a rel="nofollow" class="external text" href="https://search.worldcat.org/issn/1520-0477">1520-0477</a>. <a href="S2CID_(identifier)" class="mw-redirect" title="S2CID (identifier)">S2CID</a> <a rel="nofollow" class="external text" href="https://api.semanticscholar.org/CorpusID:124135431">124135431</a>.</cite>

</li>

<li id="cite_note-2"><a href="#cite_ref-2">^</a> <cite id="CITEREFStensrud,_David_J.2007" class="citation book cs1">Stensrud, David J. (2007). <a rel="nofollow" class="external text" href="https://books.google.com/books?id=lMXSpRwKNO8C&pg=PA56">Parameterization schemes: keys to understanding numerical weather prediction models</a>. Cambridge University Press. p. 56. <a href="ISBN_(identifier)" class="mw-redirect" title="ISBN (identifier)">ISBN</a> <bdi>978-0-521-86540-1</bdi>.</cite>

</li>

<li id="cite_note-3"><a href="#cite_ref-3">^</a> <cite id="CITEREFNational_Climatic_Data_Center2008" class="citation web cs1"><a href="National_Climatic_Data_Center" title="National Climatic Data Center">National Climatic Data Center</a> (2008-08-20). <a rel="nofollow" class="external text" href="http://www.ncdc.noaa.gov/oa/climate/conversion/swometardecoder.html">"Key to METAR Surface Weather Observations"</a>. <a href="National_Oceanic_and_Atmospheric_Administration" title="National Oceanic and Atmospheric Administration">National Oceanic and Atmospheric Administration</a>. Retrieved 2011-02-11.</cite>

</li>

<li id="cite_note-4"><a href="#cite_ref-4">^</a> <cite class="citation web cs1"><a rel="nofollow" class="external text" href="https://web.archive.org/web/20071230100059/http://weather.unisys.com/wxp/Appendices/Formats/SYNOP.html">"SYNOP Data Format (FM-12): Surface Synoptic Observations"</a>. <a href="UNISYS" class="mw-redirect" title="UNISYS">UNISYS</a>. 2008-05-25. Archived from <a rel="nofollow" class="external text" href="http://weather.unisys.com/wxp/Appendices/Formats/SYNOP.html">the original</a> on 2007-12-30.</cite>

</li>

<li id="cite_note-Krishnamurti_Annu_Rev_FM-5"><a href="#cite_ref-Krishnamurti_Annu_Rev_FM_5-0">^</a> <cite id="CITEREFKrishnamurti1995" class="citation journal cs1">Krishnamurti, T N (1995). "Numerical Weather Prediction". Annual Review of Fluid Mechanics. 27: 195–225. <a href="Bibcode_(identifier)" class="mw-redirect" title="Bibcode (identifier)">Bibcode</a>:<a rel="nofollow" class="external text" href="https://ui.adsabs.harvard.edu/abs/1995AnRFM..27..195K">1995AnRFM..27..195K</a>. <a href="Doi_(identifier)" class="mw-redirect" title="Doi (identifier)">doi</a>:<a rel="nofollow" class="external text" href="https://doi.org/10.1146%2Fannurev.fl.27.010195.001211">10.1146/annurev.fl.27.010195.001211</a>. <a href="S2CID_(identifier)" class="mw-redirect" title="S2CID (identifier)">S2CID</a> <a rel="nofollow" class="external text" href="https://api.semanticscholar.org/CorpusID:122230747">122230747</a>.</cite>

</li>

<li id="cite_note-6"><a href="#cite_ref-6">^</a> <cite id="CITEREFChaudhariLeeOh2007" class="citation book cs1">Chaudhari, H. S.; Lee, K. M.; Oh, J. H. (2007). <a rel="nofollow" class="external text" href="https://books.google.com/books?id=BQ_7vh5SrHQC&pg=PA223">"Weather prediction and computational aspects of icosahedral-hexagonal gridpoint model GME"</a>. In Kwon, Jang-Hyuk; Periaux, Jacques; Fox, Pat; Satofuka, N.; Ecer, A. (eds.). Parallel computational fluid dynamics: parallel computings and its applications : proceedings of the Parallel CFD 2006 Conference, Busan city, Korea (May 15–18, 2006). Elsevier. pp. 223–30. <a href="ISBN_(identifier)" class="mw-redirect" title="ISBN (identifier)">ISBN</a> <bdi>978-0-444-53035-6</bdi>. Retrieved 2011-01-06.</cite>

</li>

<li id="cite_note-7"><a href="#cite_ref-7">^</a> <cite class="citation web cs1"><a rel="nofollow" class="external text" href="https://web.archive.org/web/20070814044336/http://www.mmm.ucar.edu/wrf/WG4/wrfvar/wrfvar-tutorial.htm">"The WRF Variational Data Assimilation System (WRF-Var)"</a>. <a href="University_Corporation_for_Atmospheric_Research" title="University Corporation for Atmospheric Research">University Corporation for Atmospheric Research</a>. 2007-08-14. Archived from <a rel="nofollow" class="external text" href="http://www.mmm.ucar.edu/wrf/WG4/wrfvar/wrfvar-tutorial.htm">the original</a> on 2007-08-14.</cite>

</li>

<li id="cite_note-8"><a href="#cite_ref-8">^</a> <cite id="CITEREFGaffen2007" class="citation web cs1">Gaffen, Dian J. (2007-06-07). <a rel="nofollow" class="external text" href="https://web.archive.org/web/20070607142822/http://www.aero.jussieu.fr/~sparc/News12/Radiosondes.html">"Radiosonde Observations and Their Use in SPARC-Related Investigations"</a>. Archived from <a rel="nofollow" class="external text" href="http://www.aero.jussieu.fr/~sparc/News12/Radiosondes.html">the original</a> on 2007-06-07.</cite>

</li>

<li id="cite_note-9"><a href="#cite_ref-9">^</a> <cite id="CITEREFBallishKumar2008" class="citation journal cs1">Ballish, Bradley A; Kumar, V. Krishna (2008). <a rel="nofollow" class="external text" href="https://doi.org/10.1175%2F2008BAMS2332.1">"Systematic Differences in Aircraft and Radiosonde Temperatures"</a>. Bulletin of the American Meteorological Society. 89 (11): 1689. <a href="Bibcode_(identifier)" class="mw-redirect" title="Bibcode (identifier)">Bibcode</a>:<a rel="nofollow" class="external text" href="https://ui.adsabs.harvard.edu/abs/2008BAMS...89.1689B">2008BAMS...89.1689B</a>. <a href="Doi_(identifier)" class="mw-redirect" title="Doi (identifier)">doi</a>:<a rel="nofollow" class="external text" href="https://doi.org/10.1175%2F2008BAMS2332.1">10.1175/2008BAMS2332.1</a>.</cite>

</li>

<li id="cite_note-10"><a href="#cite_ref-10">^</a> <cite id="CITEREFNational_Data_Buoy_Center2009" class="citation web cs1">National Data Buoy Center (2009-01-28). <a rel="nofollow" class="external text" href="http://www.vos.noaa.gov/vos_scheme.shtml">"The WMO Voluntary Observing Ships (VOS) Scheme"</a>. <a href="National_Oceanic_and_Atmospheric_Administration" title="National Oceanic and Atmospheric Administration">National Oceanic and Atmospheric Administration</a>. Retrieved 2011-02-15.</cite>

</li>

<li id="cite_note-Hurricane_Hunters-11"><a href="#cite_ref-Hurricane_Hunters_11-0">^</a> <cite id="CITEREF403rd_Wing2011" class="citation web cs1">403rd Wing (2011). <a rel="nofollow" class="external text" href="http://www.hurricanehunters.com">"The Hurricane Hunters"</a>. <a href="Hurricane_Hunters" class="mw-redirect" title="Hurricane Hunters">53rd Weather Reconnaissance Squadron</a>. Retrieved 2006-03-30.</cite><code class="cs1-code">{{cite web}}</code>: CS1 maint: numeric names: authors list (link)

</li>

<li id="cite_note-SunHerald-12"><a href="#cite_ref-SunHerald_12-0">^</a> <cite id="CITEREFLee,_Christopher2007" class="citation news cs1">Lee, Christopher (2007-10-08). <a rel="nofollow" class="external text" href="https://www.washingtonpost.com/wp-dyn/content/article/2007/10/07/AR2007100700971_pf.html">"Drone, Sensors May Open Path Into Eye of Storm"</a>. The Washington Post. Retrieved 2008-02-22.</cite>

</li>

<li id="cite_note-13"><a href="#cite_ref-13">^</a> <cite id="CITEREFNational_Oceanic_and_Atmospheric_Administration2010" class="citation web cs1"><a href="National_Oceanic_and_Atmospheric_Administration" title="National Oceanic and Atmospheric Administration">National Oceanic and Atmospheric Administration</a> (2010-11-12). <a rel="nofollow" class="external text" href="http://www.noaanews.noaa.gov/stories2010/20100112_plane.html">"NOAA Dispatches High-Tech Research Plane to Improve Winter Storm Forecasts"</a>. Retrieved 2010-12-22.</cite>

</li>

<li id="cite_note-14"><a href="#cite_ref-14">^</a> <cite id="CITEREFStensrud,_David_J.2007" class="citation book cs1">Stensrud, David J. (2007). <a rel="nofollow" class="external text" href="https://books.google.com/books?id=lMXSpRwKNO8C&pg=PA137">Parameterization schemes: keys to understanding numerical weather prediction models</a>. <a href="Cambridge_University_Press" title="Cambridge University Press">Cambridge University Press</a>. p. 137. <a href="ISBN_(identifier)" class="mw-redirect" title="ISBN (identifier)">ISBN</a> <bdi>978-0-521-86540-1</bdi>.</cite>

</li>

<li id="cite_note-15"><a href="#cite_ref-15">^</a> <cite id="CITEREFHoughton,_John_Theodore1985" class="citation book cs1">Houghton, John Theodore (1985). <a rel="nofollow" class="external text" href="https://books.google.com/books?id=SV04AAAAIAAJ&pg=PA38">The Global Climate</a>. Cambridge University Press archive. pp. 49–50. <a href="ISBN_(identifier)" class="mw-redirect" title="ISBN (identifier)">ISBN</a> <bdi>978-0-521-31256-1</bdi>.</cite>

</li>

<li id="cite_note-Lynch_JCP-16">^ <a href="#cite_ref-Lynch_JCP_16-0">a</a> <a href="#cite_ref-Lynch_JCP_16-1">b</a> <cite id="CITEREFLynch2008" class="citation journal cs1">Lynch, Peter (2008). <a rel="nofollow" class="external text" href="https://zenodo.org/record/1326528">"The origins of computer weather prediction and climate modeling"</a>. Journal of Computational Physics. 227 (7): 3431–3444. <a href="Bibcode_(identifier)" class="mw-redirect" title="Bibcode (identifier)">Bibcode</a>:<a rel="nofollow" class="external text" href="https://ui.adsabs.harvard.edu/abs/2008JCoPh.227.3431L">2008JCoPh.227.3431L</a>. <a href="Doi_(identifier)" class="mw-redirect" title="Doi (identifier)">doi</a>:<a rel="nofollow" class="external text" href="https://doi.org/10.1016%2Fj.jcp.2007.02.034">10.1016/j.jcp.2007.02.034</a>.</cite>

</li>

<li id="cite_note-Lynch_Ch1-17"><a href="#cite_ref-Lynch_Ch1_17-0">^</a> <cite id="CITEREFLynch2006" class="citation book cs1">Lynch, Peter (2006). "Weather Prediction by Numerical Process". The Emergence of Numerical Weather Prediction. <a href="Cambridge_University_Press" title="Cambridge University Press">Cambridge University Press</a>. pp. 1–27. <a href="ISBN_(identifier)" class="mw-redirect" title="ISBN (identifier)">ISBN</a> <bdi>978-0-521-85729-1</bdi>.</cite>

</li>

<li id="cite_note-18"><a href="#cite_ref-18">^</a> Todling, Ricardo, and Stephen E. Cohn. <a rel="nofollow" class="external text" href="https://journals.ametsoc.org/view/journals/mwre/122/11/1520-0493_1994_122_2530_ssfada_2_0_co_2.xml">"Suboptimal schemes for atmospheric data assimilation based on the Kalman filter."</a> Monthly Weather Review 122, no. 11 (1994): 2530-2557.

</li>

<li id="cite_note-19"><a href="#cite_ref-19">^</a> <cite class="citation web cs1"><a rel="nofollow" class="external text" href="http://ams.confex.com/ams/91Annual/webprogram/Paper181664.html">"Abstract: Mesoscale ensemble 4DVAR and its comparison with EnKF and 4DVAR (91st American Meteorological Society Annual Meeting)"</a>. 27 January 2011.</cite>

</li>

<li id="cite_note-20"><a href="#cite_ref-20">^</a> <cite id="CITEREFYangKim2021" class="citation journal cs1">Yang, Eun-Gyeong; Kim, Hyun Mee (February 2021). <a rel="nofollow" class="external text" href="http://hfip.psu.edu/EDA2010/MZhang.pdf">"A comparison of variational, ensemble-based, and hybrid data assimilation methods over East Asia for two one-month periods"</a> (PDF). Atmospheric Research. 249: 105257. <a href="Bibcode_(identifier)" class="mw-redirect" title="Bibcode (identifier)">Bibcode</a>:<a rel="nofollow" class="external text" href="https://ui.adsabs.harvard.edu/abs/2021AtmRe.24905257Y">2021AtmRe.24905257Y</a>. <a href="Doi_(identifier)" class="mw-redirect" title="Doi (identifier)">doi</a>:<a rel="nofollow" class="external text" href="https://doi.org/10.1016%2Fj.atmosres.2020.105257">10.1016/j.atmosres.2020.105257</a>. <a href="S2CID_(identifier)" class="mw-redirect" title="S2CID (identifier)">S2CID</a> <a rel="nofollow" class="external text" href="https://api.semanticscholar.org/CorpusID:224864029">224864029</a>. Retrieved 9 November 2022.</cite>

</li>

<li id="cite_note-21"><a href="#cite_ref-21">^</a> <cite id="CITEREFBarkerLorencClayton2011" class="citation web cs1">Barker, Dale; Lorenc, Andrew; Clayton, Adam (September 2011). <a rel="nofollow" class="external text" href="https://www.ecmwf.int/sites/default/files/elibrary/2011/14950-hybrid-variationalensemble-data-assimilation.pdf">"Hybrid Variational/Ensemble Data Assimilation"</a> (PDF).</cite>

</li>

<li id="cite_note-22"><a href="#cite_ref-22">^</a> <cite class="citation web cs1"><a rel="nofollow" class="external text" href="http://www.metoffice.gov.uk/research/modelling-systems/unified-model/weather-forecasting">"Numerical weather prediction models"</a>.</cite>

</li>

<li id="cite_note-ReSeDA_report-23">^ <a href="#cite_ref-ReSeDA_report_23-0">a</a> <a href="#cite_ref-ReSeDA_report_23-1">b</a> <cite id="CITEREFBaret2000" class="citation web cs1">Baret, Frederic (June 2000). <a rel="nofollow" class="external text" href="http://w3.avignon.inra.fr/reseda/base/documents/reseda-report/00reseda-report.pdf">"ReSeDA: Assimilation of Multi-Sensor & Multi-Temporal Remote Sensing Data to Monitor Soil & Vegetation Functioning"</a> (PDF) (final report, European contract number ENV4CT960326). Avignon: <a href="Institut_national_de_la_recherche_agronomique" title="Institut national de la recherche agronomique">Institut national de la recherche agronomique</a>. p. 59. Retrieved 8 July 2019.</cite>

</li>

<li id="cite_note-24"><a href="#cite_ref-24">^</a> <cite id="CITEREFAndréGoutorbePerrier1986" class="citation journal cs1">André, Jean-Claude; Goutorbe, Jean-Paul; Perrier, Alain (1986). <a rel="nofollow" class="external text" href="https://doi.org/10.1175%2F1520-0477%281986%29067%3C0138%3AHAHAEF%3E2.0.CO%3B2">"HAPEX—MOBLIHY: A Hydrologic Atmospheric Experiment for the Study of Water Budget and Evaporation Flux at the Climatic Scale"</a>. Bulletin of the American Meteorological Society. 67 (2): 138. <a href="Bibcode_(identifier)" class="mw-redirect" title="Bibcode (identifier)">Bibcode</a>:<a rel="nofollow" class="external text" href="https://ui.adsabs.harvard.edu/abs/1986BAMS...67..138A">1986BAMS...67..138A</a>. <a href="Doi_(identifier)" class="mw-redirect" title="Doi (identifier)">doi</a>:<a rel="nofollow" class="external text" href="https://doi.org/10.1175%2F1520-0477%281986%29067%3C0138%3AHAHAEF%3E2.0.CO%3B2">10.1175/1520-0477(1986)067<0138:HAHAEF>2.0.CO;2</a>.</cite>

</li>

<li id="cite_note-25"><a href="#cite_ref-25">^</a> <cite id="CITEREFGoutorbeLebelDolmanGash1997" class="citation journal cs1 cs1-prop-long-vol">Goutorbe, J.P; Lebel, T; Dolman, A.J; Gash, J.H.C; Kabat, P; Kerr, Y.H; Monteny, B; Prince, S.D; Stricker, J.N.M; Tinga, A; Wallace, J.S (1997). "An overview of HAPEX-Sahel: A study in climate and desertification". Journal of Hydrology. 188–189: 4–17. <a href="Bibcode_(identifier)" class="mw-redirect" title="Bibcode (identifier)">Bibcode</a>:<a rel="nofollow" class="external text" href="https://ui.adsabs.harvard.edu/abs/1997JHyd..188....4G">1997JHyd..188....4G</a>. <a href="Doi_(identifier)" class="mw-redirect" title="Doi (identifier)">doi</a>:<a rel="nofollow" class="external text" href="https://doi.org/10.1016%2FS0022-1694%2896%2903308-2">10.1016/S0022-1694(96)03308-2</a>.</cite>

</li>

<li id="cite_note-26"><a href="#cite_ref-26">^</a> <cite id="CITEREFPrevotBaretChanzyOlioso1998" class="citation journal cs1">Prevot L, Baret F, Chanzy A, Olioso A, Wigneron JP, Autret H, Baudin F, Bessemoulin P, Bethenod O, Blamont D, Blavoux B, Bonnefond JM, Boubkraoui S, Bouman BA, Braud I, Bruguier N, Calvet JC, Caselles V, Chauki H, Clevers JG, Coll C, Company A, Courault D, Dedieu G, Degenne P, Delecolle R, Denis H, Desprats JF, Ducros Y, Dyer D, Fies JC, Fischer A, Francois C, Gaudu JC, Gonzalez E, Goujet R, Gu XF, Guerif M, Hanocq JF, Hautecoeur O, Haverkamp R, Hobbs S, Jacob F, Jeansoulin R, Jongschaap RE, Kerr Y, King C, Laborie P, Lagouarde JP, Laques AE, et al. (July 1998). <a rel="nofollow" class="external text" href="http://w3.avignon.inra.fr/reseda/base/watermed/documents/igarss98_e0403.pdf">"Assimilation of Multi-Sensor and Multi-Temporal Remote Sensing Data, to Monitor Vegetation and Soil: the Alpilles-ReSeDA project"</a> (PDF). Seattle, WA, USA: IGARSS'98, International Geoscience and Remote Sensing Symposium. Retrieved 8 July 2019.</cite> <code class="cs1-code">{{cite journal}}</code>: Cite journal requires <code class="cs1-code">|journal=</code> (help)

</li>

<li id="cite_note-27"><a href="#cite_ref-27">^</a> <cite id="CITEREFEiblMauserMoulinNoilhan1998" class="citation book cs1">Eibl, B; Mauser, W; Moulin, S; Noilhan, J; Ottle, C; Paloscia, S; Pampaloni, P; Podvin, T; Quaracino, F; Roujean, J.L; Rozier, C; Ruisi, R; Susini, C; Taconet, O; Tallet, N; Thony, J.L; Travi, Y; Van Leewen, H; Vauclin, M; Vidal-Madjar, D; Vonder, O.W (1998). "Comparison of the albedo derived from MOS-B and WIFS with NOAA-AVHRR". <a rel="nofollow" class="external text" href="https://hal.inrae.fr/hal-02770351/file/134612_20120110123553345_1.pdf">IGARSS '98. Sensing and Managing the Environment. 1998 IEEE International Geoscience and Remote Sensing. Symposium Proceedings. (Cat. No.98CH36174)</a> (PDF). pp. 2402–4. <a href="Doi_(identifier)" class="mw-redirect" title="Doi (identifier)">doi</a>:<a rel="nofollow" class="external text" href="https://doi.org/10.1109%2FIGARSS.1998.702226">10.1109/IGARSS.1998.702226</a>. <a href="ISBN_(identifier)" class="mw-redirect" title="ISBN (identifier)">ISBN</a> <bdi>978-0-7803-4403-7</bdi>. <a href="S2CID_(identifier)" class="mw-redirect" title="S2CID (identifier)">S2CID</a> <a rel="nofollow" class="external text" href="https://api.semanticscholar.org/CorpusID:55492076">55492076</a>.</cite>

</li>

<li id="cite_note-28"><a href="#cite_ref-28">^</a> <cite class="citation web cs1"><a rel="nofollow" class="external text" href="https://cordis.europa.eu/project/rcn/35579/factsheet/en">"ReSeDA"</a>. cordis.europa.eu. Retrieved 8 July 2019.</cite>

</li>

<li id="cite_note-29"><a href="#cite_ref-29">^</a> <cite id="CITEREFOliosoPrevotBaretChanzy1998" class="citation conference cs1">Olioso A, Prevot L, Baret F, Chanzy A, Braud I, Autret H, Baudin F, Bessemoulin P, Bethenod O, Blamont D, Blavoux B, Bonnefond JM, Boubkraoui S, Bouman BA, Bruguier N, Calvet JC, Caselles V, Chauki H, Clevers JW, Coll C, Company A, Courault D, Dedieu G, Degenne P, Delecolle R, Denis H, Desprats JF, Ducros Y, Dyer D, Fies JC, Fischer A, Francois C, Gaudu JC, Gonzalez E, Gouget R, Gu XF, Guerif M, Hanocq JF, Hautecoeur O, Haverkamp R, Hobbs S, Jacob F, Jeansoulin R, Jongschaap RE, Kerr Y, King C, Laborie P, Lagouarde JP, Laques AE, Larcena D, Laurent G, Laurent JP, Leroy M, McAneney J, Macelloni G, Moulin S, Noilhan J, Ottle C, Paloscia S, Pampaloni P, Podvin T, Quaracino F, Roujean JL, Rozier C, Ruisi R, Susini C, Taconet O, Tallet N, Thony JL, Travi Y, van Leewen H, Vauclin M, Vidal-Madjar D, Vonder OW, Weiss M, Wigneron JP (19–21 March 1998). D. Marceau (ed.). <a rel="nofollow" class="external text" href="http://w3.avignon.inra.fr/reseda/base/watermed/documents/can1.pdf">Spatial Aspects in the Alpilles-ReSeDA Project</a> (PDF). International Workshop on Scaling and Modelling in Forestry: Applications in Remote Sensing and GIS. University of Montreal, Montréal, Québec, Canada. pp. 93–102. Retrieved 8 July 2019.</cite>

</li>

<li id="cite_note-30"><a href="#cite_ref-30">^</a> <cite id="CITEREFChenWeiHsu2023" class="citation journal cs1">Chen, Shang-Ying; Wei, Jian-Yu; Hsu, Kuo-Chin (2023-10-01). <a rel="nofollow" class="external text" href="https://doi.org/10.1007/s00366-023-01897-6">"Data assimilation for real-time subsurface flow modeling with dynamically adaptive meshless node adjustments"</a>. Engineering with Computers. 40 (3): 1893–1925. <a href="Doi_(identifier)" class="mw-redirect" title="Doi (identifier)">doi</a>:<a rel="nofollow" class="external text" href="https://doi.org/10.1007%2Fs00366-023-01897-6">10.1007/s00366-023-01897-6</a>. <a href="ISSN_(identifier)" class="mw-redirect" title="ISSN (identifier)">ISSN</a> <a rel="nofollow" class="external text" href="https://search.worldcat.org/issn/1435-5663">1435-5663</a>.</cite>

</li>

<li id="cite_note-31"><a href="#cite_ref-31">^</a> <cite id="CITEREFCardenas2019" class="citation journal cs1">Cardenas, IC (2019). "On the use of Bayesian networks as a meta-modelling approach to analyse uncertainties in slope stability analysis". Georisk: Assessment and Management of Risk for Engineered Systems and Geohazards. 13 (1): 53–65. <a href="Bibcode_(identifier)" class="mw-redirect" title="Bibcode (identifier)">Bibcode</a>:<a rel="nofollow" class="external text" href="https://ui.adsabs.harvard.edu/abs/2019GAMRE..13...53C">2019GAMRE..13...53C</a>. <a href="Doi_(identifier)" class="mw-redirect" title="Doi (identifier)">doi</a>:<a rel="nofollow" class="external text" href="https://doi.org/10.1080%2F17499518.2018.1498524">10.1080/17499518.2018.1498524</a>. <a href="S2CID_(identifier)" class="mw-redirect" title="S2CID (identifier)">S2CID</a> <a rel="nofollow" class="external text" href="https://api.semanticscholar.org/CorpusID:216590427">216590427</a>.</cite>

</li>

<li id="cite_note-montabone-32"><a href="#cite_ref-montabone_32-0">^</a> <cite class="citation web cs1"><a rel="nofollow" class="external text" href="https://web.archive.org/web/20110928051446/http://www.atm.ox.ac.uk/group/gpfd/research.html#marsgcm">"Oxford Physics: Atmospheric, Oceanic and Planetary Physics: SRC: Research"</a>. July 2019. Archived from <a rel="nofollow" class="external text" href="http://www.atm.ox.ac.uk/group/gpfd/research.html#marsgcm">the original</a> on 2011-09-28. Retrieved 2011-08-19.</cite>

</li>

<li id="cite_note-hoffman-33"><a href="#cite_ref-hoffman_33-0">^</a> <cite class="citation web cs1"><a rel="nofollow" class="external text" href="http://www.eps.jhu.edu/~mjhoffman/pages/research.html">"Matthew J. Hoffman - Research"</a>.</cite>

</li>

<li id="cite_note-wrf-34"><a href="#cite_ref-wrf_34-0">^</a> <cite class="citation web cs1"><a rel="nofollow" class="external text" href="http://www.marsclimatecenter.com">"marsclimatecenter.com"</a>. marsclimatecenter.com. Retrieved 2022-04-19.</cite>

</li>

</li>

<li id="cite_note-36"><a href="#cite_ref-36">^</a> <cite id="CITEREFVetra-CarvalhoP._J._van_LeeuwenL._NergerA._Barth2018" class="citation journal cs1">Vetra-Carvalho, Sanita; P. J. van Leeuwen; L. Nerger; A. Barth; A.M. Umer; P. Brasseur; P. Kirchgessner; J-M. Beckers (2018). <a rel="nofollow" class="external text" href="https://doi.org/10.1080%2F16000870.2018.1445364">"State-of-the-art stochastic data assimilation methods for high-dimensional non-Gaussian problems"</a>. Tellus A. 70 (1): 1445364. <a href="Bibcode_(identifier)" class="mw-redirect" title="Bibcode (identifier)">Bibcode</a>:<a rel="nofollow" class="external text" href="https://ui.adsabs.harvard.edu/abs/2018TellA..7045364V">2018TellA..7045364V</a>. <a href="Doi_(identifier)" class="mw-redirect" title="Doi (identifier)">doi</a>:<a rel="nofollow" class="external text" href="https://doi.org/10.1080%2F16000870.2018.1445364">10.1080/16000870.2018.1445364</a>. <a href="Hdl_(identifier)" class="mw-redirect" title="Hdl (identifier)">hdl</a>:<a rel="nofollow" class="external text" href="https://hdl.handle.net/10754%2F630565">10754/630565</a>.</cite>

</li>

</ol></div>

<div class="mw-heading mw-heading2"><h2 id="Further_reading">Further reading</h2></div>

<ul><li><cite id="CITEREFDaley1991" class="citation book cs1"><a href="Roger_Daley" title="Roger Daley">Daley, R.</a> (1991). Atmospheric Data Analysis. Cambridge University Press. <a href="ISBN_(identifier)" class="mw-redirect" title="ISBN (identifier)">ISBN</a> <bdi>978-0-521-38215-1</bdi>.</cite></li>

<li><cite class="citation web cs1"><a rel="nofollow" class="external text" href="http://www.mmm.ucar.edu/mm5/">"MM5 community model homepage"</a>.</cite></li>

<li><cite class="citation web cs1"><a rel="nofollow" class="external text" href="https://software.ecmwf.int/wiki/display/OPTR/Data+Assimilation+Lecture+Notes">"ECMWF Data Assimilation Lecture notes"</a>.</cite></li>

<li><cite id="CITEREFIdeCourtierGhilLorenc1997" class="citation journal cs1">Ide, Kayo; Courtier, Philippe; <a href="Michael_Ghil" title="Michael Ghil">Ghil, Michael</a>; Lorenc, Andrew C (1997). <a rel="nofollow" class="external text" href="https://doi.org/10.2151%2Fjmsj1965.75.1B_181">"Unified Notation for Data Assimilation : Operational, Sequential and Variational (gtSpecial IssueltData Assimilation in Meteology and Oceanography: Theory and Practice)"</a>. Journal of the Meteorological Society of Japan. Ser. II. 75 (1B): 181–9. <a href="Bibcode_(identifier)" class="mw-redirect" title="Bibcode (identifier)">Bibcode</a>:<a rel="nofollow" class="external text" href="https://ui.adsabs.harvard.edu/abs/1997JMeSJ..75B.181I">1997JMeSJ..75B.181I</a>. <a href="Doi_(identifier)" class="mw-redirect" title="Doi (identifier)">doi</a>:<a rel="nofollow" class="external text" href="https://doi.org/10.2151%2Fjmsj1965.75.1B_181">10.2151/jmsj1965.75.1B_181</a>.</cite></li>

<li><cite class="citation web cs1"><a rel="nofollow" class="external text" href="http://www.meted.ucar.edu/nwp/pcu1/ic6/frameset.htm">"Understanding Data Assimilation"</a>. COMET module.</cite></li>

<li><cite id="CITEREFEvensen2009" class="citation book cs1">Evensen, Geir (2009). Data Assimilation. The Ensemble Kalman Filter (Second ed.). Springer. <a href="ISBN_(identifier)" class="mw-redirect" title="ISBN (identifier)">ISBN</a> <bdi>978-3-642-03710-8</bdi>.</cite></li>

<li><cite id="CITEREFLewisLakshmivarahanDhall2006" class="citation encyclopaedia cs1">Lewis, John M.; Lakshmivarahan, S.; Dhall, Sudarshan (2006). <a rel="nofollow" class="external text" href="https://archive.org/details/dynamicdataassim0000lewi">"Dynamic Data Assimilation : A Least Squares Approach"</a>. Encyclopedia of Mathematics and its Applications. Vol. 104. Cambridge University Press. <a href="ISBN_(identifier)" class="mw-redirect" title="ISBN (identifier)">ISBN</a> <bdi>978-0-521-85155-8</bdi>.</cite></li>

<li><cite id="CITEREFAschBocquetNodet2016" class="citation book cs1">Asch, Mark; Bocquet, Marc; Nodet, Maëlle (2016). Data Assimilation: Methods, Algorithms, and Applications. Society for Industrial and Applied Mathematics. <a href="ISBN_(identifier)" class="mw-redirect" title="ISBN (identifier)">ISBN</a> <bdi>978-1-61197-453-9</bdi>.</cite></li>

<li><cite id="CITEREFKalnay2002" class="citation book cs1"><a href="Eugenia_Kalnay" title="Eugenia Kalnay">Kalnay, Eugenia</a> (2002). Atmospheric Modeling, Data Assimilation and Predictability. Cambridge University Press. p. 364. <a href="Bibcode_(identifier)" class="mw-redirect" title="Bibcode (identifier)">Bibcode</a>:<a rel="nofollow" class="external text" href="https://ui.adsabs.harvard.edu/abs/2002amda.book.....K">2002amda.book.....K</a>. <a href="ISBN_(identifier)" class="mw-redirect" title="ISBN (identifier)">ISBN</a> <bdi>978-0-521-79179-3</bdi>.</cite> <code class="cs1-code">{{cite book}}</code>: <code class="cs1-code">|journal=</code> ignored (help)</li>

<li><cite id="CITEREFVetra-Carvalhovan_LeeuwenNergerBarth2018" class="citation journal cs1">Vetra-Carvalho, S.; van Leeuwen, P. J.; Nerger, L.; Barth, A.; Umer Altat, M.; Brasseur, P.; Kirchgessner, P.; Beckers, J-M. (2018). <a rel="nofollow" class="external text" href="https://doi.org/10.1080%2F16000870.2018.1445364">"State-of-the-art stochastic data assimilation methods for high-dimensional non-Gaussian problems"</a>. Tellus A: Dynamic Meteorology and Oceanography. 70 (1): 1445364. <a href="Bibcode_(identifier)" class="mw-redirect" title="Bibcode (identifier)">Bibcode</a>:<a rel="nofollow" class="external text" href="https://ui.adsabs.harvard.edu/abs/2018TellA..7045364V">2018TellA..7045364V</a>. <a href="Doi_(identifier)" class="mw-redirect" title="Doi (identifier)">doi</a>:<a rel="nofollow" class="external text" href="https://doi.org/10.1080%2F16000870.2018.1445364">10.1080/16000870.2018.1445364</a>. <a href="Hdl_(identifier)" class="mw-redirect" title="Hdl (identifier)">hdl</a>:<a rel="nofollow" class="external text" href="https://hdl.handle.net/10754%2F630565">10754/630565</a>.</cite></li></ul>

<div class="mw-heading mw-heading2"><h2 id="External_links">External links</h2></div>

Examples of how variational assimilation is implemented weather forecasting at:

<ul><li><cite class="citation book cs1"><a rel="nofollow" class="external text" href="https://www.ecmwf.int/en/elibrary/9231-part-ii-data-assimilation">Data Assimilation</a>. IFS Documentation. ECMWF. 2010.</cite> <code class="cs1-code">{{cite book}}</code>: <code class="cs1-code">|work=</code> ignored (help)</li>

<li><cite class="citation web cs1"><a rel="nofollow" class="external text" href="https://www.metoffice.gov.uk/research/weather/data-assimilation">"Data Assimilation"</a>. <a href="Met_Office" title="Met Office">Met Office</a>.</cite></li></ul>

Other examples of assimilation:

<ul><li><a rel="nofollow" class="external text" href="https://web.archive.org/web/20060717223237/http://www.cdacentral.info/">CDACentral (an example analysis from Chemical Data Assimilation)</a></li>

<li><a rel="nofollow" class="external text" href="https://web.archive.org/web/20070223072036/http://www.pdfcentral.info/">PDFCentral (using PDFs to examine biases and representativeness)</a></li>

<li><a rel="nofollow" class="external text" href="http://www.openda.org/joomla/index.php">OpenDA – Open Source Data Assimilation package</a> <a rel="nofollow" class="external text" href="https://web.archive.org/web/20110112140216/http://www.openda.org/joomla/index.php">Archived</a> 2011-01-12 at the <a href="Wayback_Machine" title="Wayback Machine">Wayback Machine</a></li>

<li><a rel="nofollow" class="external text" href="http://pdaf.awi.de">PDAF – open-source Parallel Data Assimilation Framework</a></li>

<li><a rel="nofollow" class="external text" href="http://www.data-assimilation.net/">SANGOMA New Data Assimilation techniques</a></li></ul></div><div><div class="zim-footer">

This article is issued from <a class="external text" title="Last edited on 2025-05-25" href="https://en.wikipedia.org/wiki/?title=Data_assimilation&oldid=1292186582">Wikipedia</a>. The text is available under <a class="external text" href="https://creativecommons.org/licenses/by-sa/4.0/deed.en">Creative Commons Attribution-Share Alike 4.0</a> unless otherwise noted. Additional terms may apply for the media files.

</div>

</div></div>

</div>

</main>

</div>

</body></html>